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ABSTRACT 

We present BVRIZ photometric observations of HD 209458 during the transit by its planetary com- 
panion on UT 1999 November 15 with the University of Hawaii 0.6m and 2.2m telescopes and the High 
Altitude Observatory STARE telescope. The detailed shape of the transit curve is predicted to vary 
with color due primarily to the color-dependent limb-darkening of the star, but potentially due as well 
to the effect of color-dependent opacity in the planetary atmosphere. We model the light curves and 
present refined values for the transit timing and orbital period, useful for planning future observations 
of the planetary transit. We also derive significantly improved measurements of the planetary radius, 
Rp = 1.55 ± 0.10 i?jup, stellar radius, Rs = 1.27 ± 0.05 Rq, and orbital inclination, i = 85.9 ± 0.5°. The 
derived planetary radius favors evolutionary models in which the planet has a low albedo. 
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1. INTRODUCTION 

The discovery of transits of the Sun-like star HD 209458 
by an orbiting, low-mass companion has definitively estab- 
lished the existence of extrasolar gas-giant planets (Char- 
bonneau et al. 2000; Henry et al. 2000). The identification 
of a transiting extrasolar planet system allows, for the first 
time, exploration of the physical characteristics of a planet 
outside our solar system. 

The most interesting quantity that we can derive from 
the photometric observations of the transit is the plane- 
tary radius, Rp. Although HD 209458b is less massive than 
Jupiter, its radius is significantly larger, a result of slowing 
of the planetary contraction due to exposure to high stel- 
lar insolation early after the planet's formation (Burrows 
et al. 2000). The precise value of Rp is dependent upon 
the planetary Bond albedo. A, a fundamental parameter 
in the thermodynamics of the planetary atmosphere. A 
precise measurement of Rp may allow an estimate of A 
(Burrows et al. 2000), though the published uncertainty 
in Rp (Mazeh et al. 2000) is too great for a meaningful 
constraint. 

The photometric observations of the transit published 
to date have been in the R band (two full transits; Char- 
bonneau et al. (2000)) and in a Stromgren b + y band 
(one half transit; Henry et al. (2000)). In this Letter, 
we present observations in five wavelength bands of the 
planetary transit on UT 1999 November 15. As we de- 
scribe below, we have observed color-dependent variations 
in the transit shape due to the stellar limb-darkening. We 
present a likelihood analysis of the data that improves the 
estimates of the planetary radius and orbital inclination. 



More generally, we outline a method that may be used 
to analyze future high-precision, multicolor data and ex- 
ploit the color-dependence in order to break the degener- 
acy shared between the planetary and stellar parameters. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. University of Hawaii 2.2m Data 

Our multicolor observations of the HD 209458 planetary 
transit were made with the University of Hawaii 2.2m tele- 
scope on Mauna Kea under photometric skies. The data 
were taken with the Tek 2048^ CCD at 0.22" per pixel. 
The brightness of HD 209458, V = 7M,B-V = 0.58 
(H0g et al. 2000)^, forced us to defocus the telescope in 
order not to saturate the detector. As there were no stars 
of comparable brightness in the field of view, we took im- 
ages of a nearby comparison star before and after each set 
of HD 209458 observations. Each observation sequence 
on target consisted of five 1 s exposures in each of John- 
son V and Kron-Cousins R and /, and five 5 s expo- 
sures with a Z filter. We also observed comparison stars 
that bracketed HD 209458 in airmass, alternating between 
HD 210483 {V = 7.58, B~V = 0.61) and HD 208156 
{V = 8.11,-8 — F — 0.50). The comparison star obser- 
vations were taken in the same manner as those of the 
target star. 

We performed aperture photometry on the bias- 
corrected and flat-fielded CCD images, using a 60" digital 
aperture. The magnitudes of the comparison star were 
plotted against airmass to define an extinction correction. 
We continued observations as far as airmass 2.1, but a lin- 
ear extinction correction was sufficient to fit the data. We 
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also corrected for color-dependent extinction by interpo- 
lating between the extinction corrections derived for the 
comparison stars (which straddled HD 209458 in color). 

We converted the derived magnitues into relative fluxes 
and each set of five consecutive observations (in each fil- 
ter) was binned to produce a single average value, with an 
uncertainty given by the error in the mean. The photom- 
etry was limited primarily by scintillation noise, and the 
achieved relative precision for the time series was 1.5 x 10~^ 
(1.6 mmag), 9 x 10^'' (1.0 nimag), 1.1 x 10^"^ (1.2 mmag), 
and 7 x lO"** (0.8 mmag), for the V, R, I, and Z data, 
respectively. For each band, the off-transit flux level was 
derived from the statistically-weighted average of the last 
three observations, and time series in each color were nor- 
malized to this level to yield the relative fluxes used in the 
modeling. Applying the reduction procedure to the com- 
parison star observations yielded constant light curves as 
expected. 

2.2. University of Hawaii 0.6m Data 

We obtained Johnson B observations of HD 209458 with 
the University of Hawaii 0.6m telescope, in combination 
with a three channel photometer, similar to the one de- 
scribed in Kleinman, Nather, & Phillips (1996). Using 
identical photomultiplier and amplifier-discriminator de- 
tectors this device allows high-speed aperture photometry 
to be simultaneously conducted on three selectable regions 
of the sky, thus enabling continuous monitoring of a target 
star, comparison star and sky background. 

HD 209458 was continuously monitored using 10 s in- 
tegrations through a Johnson B filter. A nearby (fainter) 
star was used for comparison and guiding, with all chan- 
nels (including sky) using 30" apertures. A B filter was 
not available for the comparison star channel, so its effec- 
tive filter was determined by the atmosphere to the blue 
and a dichroic beam-splitter to the red (sending the red 
light in this channel to a CCD guider), resulting in a bluer 
and broader passband than the target channel. We con- 
sequently obtained an improved atmospheric absorption 
model for the final data reduction by observing HD 209458 
off-transit on a subsequent night (UT 1999 November 17), 
during virtually identical photometric conditions over the 
same airmass range (1.1 to 1.8) and with the same instru- 
mental configuration. The HD 209458 sky-corrected in- 
strumental magnitudes were corrected for color-dependent 
extinction with a linear airmass model. These were con- 
verted to relative fluxes and binned into 4 minute time 
bins, with normalization by the weighted average of the fi- 
nal two data bins, which occurcd off-transit. The resulting 
rms relative flux residuals were 7 x 10^** (0.8 mmag). 

2.3. STARE Telescope Data 

Additional R observations of HD 209458 were made with 
the High Altitude Observatory STARE telescope (Brown 

& KoHnski 1999'^; Brown & Charbonneau 2000) on the 
same night, though the target star set at mid-transit. The 
observations and reduction of these data were identical 
to those described by Charbonneau et al. (2000), though 
the precision is significantly lower than those observations 
since the star was observed at higher airmass. The data 
were binned into 5 minute sets, and converted to relative 

'^Available at http://www.hao.ucar. edu/public/research/stare/stajre. 



fiuxes by normalizing by the statistically- weighted average 
of the off-transit data. The achieved relative flux residual 
rms increased from 3.7 x 10~^ (4 mmag) to 6.3 x 10"'^ (7 
mmag) over the course of the night as the target set. 

Our final multicolor photometry is displayed in Figure 
1, along with our best-fit model, described in the next 
section. 

3. ANALYSIS 

3.1. Improvement to Orbital Parameters 

Assuming the transit is symmetric, we find the time for 
the center of the transit is T^ = 2451497.797 ± 0.002 HJD. 
The two previous observations of the full transit (Char- 
bonneau et al. 2000) found % = 2451430.823 ±0.003 HJD 
and Tc = 2451437.873±0.003 HJD, which occured, respec- 
tively, 19 and 17 orbits earlier than the one we report here. 
From these three measurements of Tc, we calculate an or- 
bital period of P = 3.52495 ± 0.0003 days. We note also 
that this period is consistent with the more precise state- 
ments of the period ofP = 3.524736±0.000045 days (3.9 s) 
(Castellano et al. 2000) and P = 3.524739 ±0.000014 days 
(1.2 s) (Robichon & Arenou 2000) derived from Hipparcos 
archive photometry. 

3.2. Analysis of the Multicolor Transit Curve 

Derivation of parameters from transit observations in a 
single bandpass is plagued by a degeneracy: the transit 
curve is described primarily by its depth and duration, 
yet it is always possible to fit these with a larger (smaller) 
planet provided the stellar radius is also increased (de- 
creased) and the orbital inclination, i, is decreased (in- 
creased). However, by observing the relative depth and 
shapes of the transit curves in different bands, this de- 
generacy can be lifted. Measurement of the relative flux 
at a given time during the transit can be combined with 
a limb-darkening model for the star, /a(m)j to derive /i, 
the cosine of the angle between the normal to the stellar 
surface and the line of sight, thus constraining the inclina- 
tion. With this constraint, the transit duration constrains 
the stellar radius, which when combined with the transit 
depth (a measurement of the ratio Rp/ Rs) constrains the 
planetary radius. For example, if one observes that the 
transit is deeper in R than in B, one would infer that a 
larger planet passes close to the limb of a larger star, be- 
cause the star is relatively brighter in the red towards the 
limb. 

Because we can observe only the flux integrated over the 
stellar disk, we proceed as follows. R. Kurucz (personal 
communication) generated model stellar intensities as a 
function of wavelength and /i. For each passband, we mul- 
tiplied the model intensities by a response function (see be- 
low) and integrated over wavelength to produce intensity 
as a function of /i. For th(^ 2.2m obsc^rvations, the response 
function was the combination of the mirror, CCD, and fll- 
ter responses, whereas for the 0.6m and STARE data, we 
used only the filter response. We then produced, for each 
passband, C, a grid of model light curves, Mc{Rp, Rs,i), 
as a function of the planetary radius, the stellar radius, 
and the orbital inclination. The light curve model follows 
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the description presented in Sackett (1999) and Charbon- 
neau et al. (2000), taking into account the finite size of 
the planet and the limb-darkening of the star. Since we 
integrate over the whole stellar disk and a wide passband, 
the light curves are relatively insensitive to small errors in 
the Kurucz model; we have confirmed this with fits using 
parameterized limb- darkening models. 

Because all data points of the light curve have been nor- 
malized by the statistically weighted average of the off- 
transit points, there is a potential systematic uncertainty 
in the transit depth, as there are fewer off-transit points 
than on-transit points. Thus we introduced a fourth pa- 
rameter into the model, a scale factor for the data, /, 
equivalent to allowing for an unknown magnitude zero- 
point. This nuisance-parameter was subsequently inte- 
grated out (see below), and this integration serves to ac- 
count properly for the systematic uncertainty introduced 
by the paucity of out-of-transit data by appropriately 
broadening the likelihood distributions. 

We adopt a Bayesian approach for the analysis; in a 
given passband the likelihood of a particular model is 

PciRp, Rs, i, f) oc exp (^-y) (1) 
with 

2 f Dcjtn) - fMcitn;Rp,Rs,i) y 

^ h V ^C{tn) J ' ^' 

where Dc{tn) are the observed data with uncertainties 
(Jc{tn), and the p are the prior probabilities (priors; see 
below). We evaluated pc separately in each passband 
for our grid of light curve models spanning the range of 
{Rp, Rs, i, /} where the likelihood was non-negligible. 

As we wish to measure the parameters independently of 
the previous results, we assume we have no prior knowl- 
edge of Rp and i. However, we can estimate Rg as fol- 
lows. With the photometry as above, the Hipparcos par- 
allax of HD 209458 yields My = 4.27 ± 0.11 (Ferryman 
et al. 1997). Applying a bolometric correction (Bessell, 
Castelli, & Plez 1998) of —0.06, appropriate for the stellar 
type and metallicity (Mazeh et al. 2000), we can calculate 
the stellar bolometric flux. Then from the measured effec- 
tive temperature of 6000 K (Mazeh et al. 2000), we derive 
Rs = 1.18i?0. There is a 10% uncertainty in the visual 
luminosity and a 1% uncertainty in the effective temper- 
ature, and we estimate a few percent uncertainty in the 
bolometric correction. Combining these implies an 8% un- 
certainty in the stellar radius. Thus, in our analysis of the 
transit curves we adopt a Gaussian prior for the stellar 
radius, p{Rs) corresponding to Rg = 1.18 ± O.O9i?0, and 
uniform priors p{Rp), p{i) and p{f), on the other model 
parameters. Different choices for the exact shape of the 
prior on the stellar radius did not significantly affect the 
results. 

We derived pc independently for each of the six pass- 
bands: the UH G.6m B data; the UH 2.2m V, R, I and Z 
data; and the STARE R data. The minimum reduced is 
very near unity for four of the six passbands, B : xt = 0-97 
(with 39 degrees of freedom), V : xl = 1-16 (13), STARE 
R:xl^ 1-03 (49), I ■.xl = 1-04 (13). The two exceptions 
are the UH 2.2m R and Z data, which have minimum x^ 
of 2.85 and 2.97, respectively for 13 degrees of freedom 
each. As these are the two passbands with the smallest 



statistical uncertainty, we adopt the plausible hypothesis 
that the poor fits are caused by an unidentified systematic 
uncertainty (at the level of ~ 1 mmag) , rather than a fail- 
ing of the light curve models. Data of higher precision will 
resolve this issue, but for the time being we have scaled 
the uncertainties for the UH 2.2m R and Z data to give a 
minimum reduced x^ of unity, thus broadening the derived 
likelihood distributions. 

To derive combined constraints from the individual like- 
lihoods, we first integrate each pc over the nuisance pa- 
rameter /, leaving correctly broadened three-dimensional 
likelihood distributions over {Rp, Rs,i}- Multiplying these 
distributions together (but with the priors included only 
once) yields our combined likelihood; the peak of this dis- 
tribution is our best-fit model and the light curves cor- 
responding to this model are shown as the solid lines in 
Figure 1. We note that the uncertainties in our derived 
parameters are highly correlated and to underscore this 
important fact, we present contours of the two-dimensional 
likelihood distributions (integrated over the third variable) 
in Figure 2. To correctly account for these correlations 
in determining the individual likelihoods for each of the 
parameters separately, we integrate the three-dimensional 
distribution over the other two parameters. These one- 
dimensional distributions are well- fit by Gaussians, with 
Rp = 1.55 ± 0.10 iijup, Rs = 1.27 ± 0.05 i?0, and i = 
85.9 ±0.5°. 

4. DISCUSSION AND CONCLUSION 

Our derived parameters are consistent with the results 
of Mazeh et al. (2000), with our multicolor data roughly 
halving the uncertainties presented there, and favoring a 
slightly larger planet and star. Combined with the orbital 
data of Mazeh et al. (2000), our results imply a mean plan- 
etary density of p = 0.23 ± 0.05 gcm~^. 

We are also in a position to make a comparison between 
the observations and models for HD 209458b. Burrows 
et al. (2000) present evolutionary tracks of radii versus 
time for close-in extrasolar giant planets, taking into ac- 
count the effects of high stellar insolation in retarding the 
planet's contraction. Previous measurements of the radius 
of HD 209458b were not precise enough to differentiate be- 
tween models with low and moderate albedos (models A 
and B of Burrows et al. (2000) in their Figure 1, with 
Bond albedos of 0.0 and 0.5, respectively). Our more pre- 
cise measurement favors a larger planet, and is in very 
good agreement with model A, which has a radius close to 
1.5 iijup (at the approximate age of the system of several 
Gyr). Our derived radius disfavors model B, for which 
the radius is predicted to be near 1.3 i?jup. This result 
provides evidence for a low albedo for HD 209458b. 

It is also of interest to see if we can detect any differences 
in the planetary radius as a function of color. The value 
of Rp as derived from transit observations may vary with 
the passband, since the planet would appear slightly larger 
when observed at wavelengths where the atmosphere con- 
tains strong opacity sources (Brown & Charbonneau 2000; 
Burrows et al. 2000; Seager & Sasselov 2000). To explore 
this possibility, in Figure 3 we present the peak positions 
of the one-dimensional conditional likelihood distributions 
for the planetary radius in each passband, with the stellar 
radius and inclination fixed at their best-fit values. The 
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uncertainties shown then do not incorporate the correla- 
tions exhibited in Figure 2, but they are correct relative 
to each other. For instance, if the stellar radius were in- 
creased, all of the points on Figure 3 would increase to- 
gether. Even so, the data arc not precise enough to show 
evidence for any variation in the planetary radius with 
wavelength. Future observations should be able to im- 
prove the measurement precision, and provide a tool for 
learning about the planetary atmosphere. 

The discovery of an cxtrasolar planet transiting its host 
star has heralded a new era in a young field. While pre- 
cision radial velocities reveal much about orbital charac- 
teristics, the addition of transit observations yield under- 
standing about the physical characteristics of the plan- 
ets themselves. We eagerly anticipate further exploration 



of extrasolar planets with continued observations of this 
shadowy nature. 
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Fig. 1. — Multicolor observations of the HD 209458 planetary transit on UT 1999 November 15, from the University of Hawaii 2.2m (circles) 
and 0.6m (crosses) telescopes and the High Altitude Observatory STARE telescope (triangles). Uncertainties in the B data are smaller than 
the symbols. The error bars do not include the systematic uncertainties in the zero points, though these uncertainties were included in the 
fitting procedure. The solid curves are our best-fit model for the combined data set. 
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Fig. 2. — Contour plots of the marginalized two-dimensional likelihood distributions, illustrating that all three parameters are highly 
correlated. The peak position is shown by a cross, and the contours displayed correspond to 68.3, 95.4 and 99.7% confidence levels. 
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Fig. 3. — Planetary radii derived independently for each of the six data sets, with the stellar radius and orbital inclination fixed at 
Rs = 1.27-Rjup and i = 85.9°, respectively. The symbols arc as in Figure 1, with the points placed at the cfTcctivc wavelength of each 
bandpass and the horizontal error bars indicating where the transmission falls to half the peak value. The vertical error bars do not include 
the uncertainty common to all the points from the uncertainties in the stellar radius and orbital inclination. The dotted line shows the the 
best-fit planetary radius with all the data combined. 



